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We consider the algorithmic problem of selecting a set of target nodes that cause the biggest 
activation cascade in a network. In case when the activation process obeys the diminishing returns 
property, a simple hill-climbing selection mechanism has been shown to achieve a provably good 
performance. Here we study models of influence propagation that exhibit critical behavior, and 
where the property of diminishing returns does not hold. We demonstrate that in such systems, the 
structural properties of networks can play a significant role. We focus on networks with two loosely 
coupled communities, and show that the double-critical behavior of activation spreading in such 
systems has significant implications for the targeting strategies. In particular, we show that simple 
strategies that work well for homogeneous networks can be overly sub-optimal, and suggest simple 
modification for improving the performance, by taking into account the community structure. 

PACS numbers: 89.75.Hc 



I. INTRODUCTION 

Much recent research has focused on understanding 
how structural properties of networks affect their dy- 
namical behavior [H [21 El H] . For instance, it has been 
established that critical behavior of epidemic models on 
random Erdos-Renyi graphs are absent in certain scale- 
free networks [5l [6]. This fact has significant implica- 
tions for a number of important applications. In par- 
ticular, while scale-free networks are generally robust to 
random breakdowns, they can be highly vulnerable to in- 
tentional attacks that target highly connected nodes [Tj. 
This suggests that immunization strategies for such net- 
works should take into account the inherent heterogene- 
ity in the degree distribution. Indeed, it has been shown 
that targeted immunization based on the nodes connec- 
tivity hierarchy can significantly lower the networks vul- 
nerability to epidemic attacks [8J. 

Here we consider a related problem of maximizing in- 
fluence propagation in networks, by targeting certain in- 
fluential nodes that have the potential to influence many 
others. This problem has attracted some recent attention 
due to potential applications in viral marketing, which is 
based on the idea of leveraging existing social structures 
for word-of-mouth advertising of products [21 [TT] . 
From the algorithmic standpoint, this selection problem 
can be stated as follows [Hd^]: Given a social network, 
an influence model, and a set of nodes S, let a{S) be 
the expected number of nodes that will be activated by 
the end of the influence propagation process. Then, for 
a given budget n, the influence maximization problem 
is concerned with finding the set S of size n that max- 
imizes the return cr(S). While this problem is known 
to be NP hard for the many influence models, several 
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approximate methods have been developed. An impor- 
tant result established in [12] states that for a class of 
models for which the return function is sub-modular, a 
simple hill-climbing algorithm, which works by greedily 
selecting the next best candidate node, yields a solution 
which is guaranteed to be within ~ 63% of the opti- 
mal. Sub-modularity of the return function means that 
a{S U {oj}) - a{S) > a{T U {w}) - cr{T) for any node uj 
and any 5 C T. In other words, the expected return for 
targeting a node diminishes with the number of targeted 
nodes. 

While it is quite safe to assume that the diminish- 
ing returns property is satisfied in saturated, or near- 
saturated, markets, those models might fail to capture 
the dynamics of emerging markets, where the condition of 
the sub-modular growth can be violated. Indeed, many 
economical and social phenomenon are better described 
in terms of critical phase transitions, where a huge growth 
is observed only after some threshold conditions are met. 
Here we are interested in this latter case. Our main re- 
sult is that in such critical systems, the structural proper- 
ties of networks can play a significant role in the dynam- 
ics of the influence propagation. Consequently, selection 
strategies that do not account for those structural prop- 
erties might produce vastly sub-optimal results. 

To be more speciflc, let us focus on the so called linear 
threshold models (LTM) [TSl where a node is acti- 
vated whenever the fraction of its active neighbors ex- 
ceeds some pre-defined threshold, J2jeAfi — Here 
A/i is the set of active neighbors of node i, Wij is the 
normalized weight of the link between the nodes i and j , 
J2j = 1, and < < 1 is the activation threshold 
for the node i. Usually, 0j-s are assumed to be random 
variables refiecting the uncertainty about individuals. 

Consider the unfolding of the LTM dynamics on the 
hypothetical influence graphs depicted in Figure [l] In 
Figure [l] (a), the nodes in the upper row are connected in 
pairs, and each link has a weight 1/2. Each node from the 
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FIG. 1: Two hypothetical networks illustrating the limitation 
of the hill-climbing algorithm. In both networks, each node 
from the bottom (upper) row affects all the upper (lower) 
nodes with weight 1/2K {1/N). In (b), each of the K nodes 
in the lower row is influenced by its two preceding neighbors 
with weights 1/4. 



bottom (upper) row affects all the upper (lower) nodes 
with weight 1/2K {l/N). Assume fixed thresholds 1/2 
for all the nodes . If one follows the hill-climbing algo- 
rithm, then it is easy to see that only the upper nodes 
will be selected. Thus, after targeting n nodes, the num- 
ber of active nodes is exactly 2n. li n < K, then this 
is indeed an optimal solution to influence maximization 
problem. However, for n > if this solution is clearly sub- 
optimal as targeting the bottom K nodes would activate 
all the N + K nodes in the network. The sub-optimality 
is even more dramatic for the network depicted in Fig- 
ure [T](b), which is obtained from Figure [T] (a) by adding 
2K links so that each of the lower nodes is now influ- 
enced by its two proceeding neighbors with weights 1/4. 
Assuming a threshold 1/2 for all the nodes, one observes 
that the greedy selection policy will again result in a final 
active set of size 2n. However, a simple inspection shows 
that if one activates two neighboring nodes from the bot- 
tom row, then it will cause a global cascade among the 
lower nodes, which will consequently propagate to the 
upper nodes and activate them as well. This suggests 
that for large N, greedy selection mechanism produces 



vastly sub-optimal solution. 

While the examples above seem peculiar, the main 
claim of our paper is that the underlying effect is rather 
general and present in more realistic models as well. In- 
deed, the two contributing factors to the behavior de- 
scribed above are the critical nature of the activation dy- 
namics, and the structural heterogeneity of the network. 
The criticality is manifested by the fact that there is a 
threshold number so that for n < influence prop- 
agation is localized, whereas for n > ric the activation 
spreads throughout all (or almost all) the nodes in the 
network. And by the structural heterogeneity we mean 
different and heterogeneous linkage-patterns among the 
nodes. A large class of networks that fit this description 
are networks with well-defined communities |15l I16j . In 
particular, here we focus on networks that are composed 
of a relatively small, tight community that is connected 
with a larger population of nodes (see Figure l2l . 




FIG. 2: Schematic illustration of a bi-community network. 



II. ACTIVATION DYNAMICS ON 
BI-COMMUNITY ERDOS-RENYI GRAPHS 

We have previously analyzed activation dynamics on 
networks composed of two loosely coupled Erdos-Renyi 
graphs |17j . To make this paper self-contained, we be- 
low provide our analysis. Instead of using the traditional 
linear threshold model, we focus on a modified version, 
where the threshold condition is applied not to the frac- 
tion of active neighbors, but their number. Thus, a node 
is activated whenever the number of its active neighbors 
is greater or equal a predefined threshold, h. The reason 
for modifying the model is that the consequent analysis 
is simple. Furthermore, one can argue that the modified 
model might more plausible from the social choice stand- 
point: Indeed, it is hard to imagine that, while making a 
decision based on the opinion of few friends, one "weighs" 
the advice by the overall number of friends. We stress, 
however, that our main results are valid for the fractional 
threshold model as well, provided that it possesses criti- 
cal behavior. 

Let us first focus on a single Erdos-Renyi graph with 
an average connectivity z. Let Ph denote the fraction of 
nodes with threshold h, and let po be the budget, i.e., the 
fraction of targeted (initially activated) nodes. In this 
section we consider random uniform targeting, so that 
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each node has a probability po to be targeted. Following 
the same line as in [17] it can be shown that the fraction 
of activated nodes at the end of the cascading process 
satisfies the following transcendental equation: 



= 1 



{l-po)Y,PhQ{h-zp') 



(1) 



h=0 



where Q{n,x) = X]fc<n ^~^^'^/^' the regularized 
gamma function. To understand this expression, note 
that in the steady state, a node with a threshold h is in- 
active if it is connected with less than h active neighbors, 
and it is not among the initially targeted set. The former 
happens with probability Q{h;zp''), while for the latter 
this probability is {1 — po), hence yielding Equation [l] 

For a fixed connectivity z, the solution of the Equa- 
tion [l] depends on the budget po, as well as on the thresh- 
old distribution function Pf^. Let us elaborate on the 
latter dependence in more details. First of all, we as- 
sume that Po = 0, i.e., there are no nodes that activate 
spontaneously, aside from the initially targeted nodes. 
Furthermore, simple inspection shows that the dynam- 
ical properties of the model depend on the fraction of 
nodes with threshold h = 1, Pi. Following [T3], we call 
these nodes vulnerable since they will activate whenever 
one of their neighbors is active. Clearly, if the fraction 
of the vulnerable nodes is sufficiently large, a single node 
might trigger a global cascade throughout the network. 
In particular, a global cascade will happen whenever the 
vulnerable nodes form a giant connected component |14j , 
which, for the random Erdos-Renyi graphs translates 
into Piz = c > 1. Below we consider the case when 
Pi is either zero, or sufficiently small. Pi ^ 1/z, so that 
for a network of size N, the number of nodes required to 
cause a global cascade must be of order 0{N) as N ^ oo. 

For the latter case, the analysis of Equation[l]yields the 
following observation: For a given connectivity z, there is 
a critical fraction pc such that for po < pc the activation 
process is localized, while for po > pc activation spreads 
to all the nodes in the network. One can obtain the 
following expression for the critical density: 
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where xq satisfies the following equation: 
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Here Dh — 'Ylii<h '^^ the cumulative distribution func- 
tion for the activation thresholds. 

Consider now the activity spreading in two coupled 
Erdos-Renyi networks of sizes Na and iVf, as depicted in 
Figure [2j with connectivities Zaa,Zbb within the groups, 
and Zab — ZbaNb/Na across the groups. Assume that 
the cascading process in group A is not affected by cross- 
group links, so that the activation for A nodes is governed 



by the Equation [T] For the B nodes, the activation dy- 
namics is given by a similar equation, with the only dif- 
ference that it is affected by the presence of active A 
nodes, and the steady state fraction of active B nodes 
satisfies the following equation: 
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where is the steady state fraction of active A nodes. 
Thus, the presence of active A nodes facilitates the acti- 
vation of B nodes, and the effect depends on the across 
the group connectivity Zba- Specifically, if Zba is very 
small, then, in order to achieve a global activation in 
group B, one needs to target fraction of B nodes above a 
certain threshold pb^c- However, even below the thresh- 
old, there is a possibility of a global cascade in group 
B if the across the group connectivity Zba is sufficiently 
large. Indeed, a simple analysis showsjlT] that for a fixed 
within-group connectivity Zbb, there is a critical connec- 
tivity z^^ so that for Zba > z^^ the activation will propa- 
gate from group A to group B and cause a global cascade. 



A. Influence Maximization in bi— community 
Erdos-Renyi graphs 

The analysis above suggest that discarding the com- 
munity structure might result in sub-optimal solution to 
the influence maximization problem. Indeed, since the 
critical number of nodes necessary to cause a cascade for 
a given connectivity grows linearly with the network size, 
it might be more beneficial to target the smaller group 
first and cause an activation cascade in that group. Af- 
terwards, the activation will propagate through the larger 
network, provided that the density of links between the 
groups is sufficiently strong. Strictly speaking, the anal- 
ysis above applies to the random targeting strategies. 
However, one might expect a similar reasoning to hold 
for the greedy selection heuristics as well. Indeed, below 
we validate this hypothesis for synthetic random graphs, 
which which are similar to those used in the evaluation of 
community finding algorithms |18j . Namely, we assume 
that the network is composed of L groups, with Nl nodes 
in each. Each pair of nodes within and across the same 
groups are linked with probability pi„ and Pouti respec- 
tively, with corresponding connectivities Zi„ = PmNm 
and Zout — Pout{N — Nl). We assume that one of those 
L groups constitute the group A, while the remaining 
L — 1 communities form B. 

We tested the greedy selection algorithm for with both 
integer and fractional versions of the linear threshold 
model. The conventional greedy selection works as fol- 
lows: Starting from the empty set S* = 0, the algorithm 
finds a single node that causes the largest cascade, adds it 
to S (if there are many candidates nodes, then one needs 
a tie-breaking mechnaism: Here we choose the node with 
the highest connectivity). This process is then repeated 
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lected. 

The same picture holds for the fractional-threshold 
mode as shown in Figure |4](b), where we again compare 
both algorithms. The fractional thresholds are uniformly 



distributed on the interval [9min,Qmax\, with i 



0.1, 



(^max = 0.5 used here. Again, we observe that for both 
small and large values of n, both methods have a similar 
performance, while for an intermediate values of n, the 
strategy that selects A nodes are superior. 



III. INFLUENCE MAXIMIZATION ON 
SCALE-FREE GRAPHS 
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FIG. 3: Comparison of target-selection strategies for integer 
(upper panel) and fractional (lower panel) threshold models. 
We used L = 10, and Nl ~ 500, and the total network size 
is N = 5000. The connectivities were set to Zi„ = Zout = 8, 
and the thresholds were chosen randomly (uniformly) from 
the intervals shown in the inset. 



n times. We compared this simple hill-climbing scheme 
with another one, which works exactly the same way, 
but now the candidate nodes are selected only from the 
smaller community A. In the following, we differentiate 
the latter algorithm by a subscript A. 

In Figure [3] (a) we plot the fraction of activated nodes 
against the budget n, for the two selection strategies, 
and for the integer threshold model. The connectivities 
are set to Zi„ = -Zout = 8, and the thresholds were chosen 
randomly and uniformly from the interval [2, 5]. One can 
see that the strategy of targeting nodes from the smaller 
community is generally more efficient, as it achieves a 
global cascade with a significantly lower budget n. More 
precisely, for small and large values of n, both methods 
have a similar performance. However, there is a win- 
dow [n^ , rig] , within which the selection of A nodes is 
clearly superior. Recalling the analysis from the previ- 
ous section, it is clear that corresponds to the critical 
threshold for which the activation spreads throughout 
group A, and then spills into the rest of the network. If 
one targets nodes from the general population, on the 
other hand, this critical effect does not come into play 
until later, when larger number of nodes have been se- 



Real-world networks have statistical characteristics 
that significantly deviate from the Erdos-Renyi model. 
In particular, many networks exhibit power-law degree 
distribution. Below we examine influence maximization 
in such networks. We show that the sub-optimality of 
the simple greedy algorithm persists whenever the con- 
nectivity variance is limited, so that activation dynamics 
still demonstrates critical behavior. 

Let us first focus on the activation dynamics in a single 
population. Furthermore, for the sake of simplicity, let 
us assume that all the nodes have the same activation 
threshold H. We consider a network with degree distri- 
bution pk — ck~'^ , m < k < M, where m and M are 
lower and upper cutoffs, respectively, and c is a normal- 
izing constant. 



(7-1) 



(3 = m/M 



(5) 



Let pk be the fraction of nodes with connectivity k that 
are active at the end of the activation process, and let 
Po.k be the fraction of targeted A:-nodes. Using similar 
arguments as in the Erdos-Renyi case, one can show that 
PkS satisfy the following equations {k — m,m + 1, ..M): 



H-l 



Pfc = l-(l-po.fc)^P(m|fc). 



(6) 



Here P{m\k) is the probability that m out of k edges 
leaving from a vertex point to an activated node. Let 9 
be the probability that a randomly chosen edges leads to 
an active node. Then P(m\k) is a binomial distribution 
that for large k can be approximate by the Poisson dis- 
tribution with a mean 9k. Furthermore, for uncorrelated 
networks considered here, 9 can be written as 



'-Y^kp^ 



kPk, 



(7) 



where 2 = kp^ is the average connectivity. To un- 
derstand this expression, note that the probability that 
a randomly chosen edge leads to a node with degree k is 
proportional to kp\^ (for uncorrelated networks) , and the 
probability that this node will be active is simply p^. 
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Combining Equations [6] and [7] we obtain the following 
self-consistent equation for 6 in the continuos approxi- 
mation 
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To proceed further, we need to specify the targeting 
function po.fc- We considered two cases - random se- 
lection /5o,fc = Po = const; and maximum degree (MD) 
selection heuristics po,fe = ©(^^"^0)1 where @ is the step 
function, and the cutoff mg is found from the budgeting 



constraint J^^^^Pk = Poi which yields 



mo 



(9) 



Examination of Equation |8] for both strategies can be 
summarized as follows. First of all, it is easy to see that 
for sufficiently dense networks, = 1 is always a so- 
lution. Thus, for sufficiently large po the steady state 
corresponds to a fully activated network. Furthermore, 
for 7 > 3, there is a critical fraction pc below which 
another solution appears, as shown schematically in Fig- 
ure |4] In the region 2 < 7 < 3, the critical behavior is 
suppressed if there is no upper connectivity cutoff, f3 = 0. 
Namely, for any finite po, the network is fully activated 
at the end of the cascading process, for arbitrary acti- 
vation threshold H. This is due to the infinite second 
moment of the connectivity distribution. The criticality 
is recovered, however, if one introduces an upper cutoff. 
More details axe provided in the Appendix. 
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FIG. 5: Comparison of target-selection strategies for scale- 
free network with 7 = 2.5. We used m = 10, and Na = 1000, 
and the total network size is = 5000. 
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FIG. 4: Graphical illustration of the steady state equation for 
the random targeting strategy. 

Next, we examine the predictions of the above analysis 
for the influence maximization problem in networks with 
power-law degree distribution. The networks were gener- 
ated according to the configuration model [121 HO] i with a 
slight modification to account for a community structure. 



Namely, we assigned each node to one of two communi- 
ties. Then, if the generated candidate edge was linking 
nodes in different communities, that edge was rejected 
with probability Thus, = corresponds with a 
single-community scale-free network, while pi„ = 1 cor- 
responds to two completely disjoint networks. We choose 
•pin such that the network has a well-defined commu- 
nity structure, but at the same time, the number of links 
across the communities is sufficiently large so that the 
cascade can spread from one community to the other. 

We examined the impact of community structure on in- 
ffuence maximization for several targeting strategies and 
for varying connectivity cutoff M . For relatively smaller 
cutoffs (Af ~ m), the behavior of the activation spread- 
ing should be similar to the results for the random Erdos- 
Renyi graphs, while for larger M it recovers the scale free 
characteristics. In Figure [5] (a) and (b) we compare the 
two variants of the maximum degree heuristics and the 
greedy strategy, respectively, by plotting the minimum 
expected budget one needs to achieve a global cascade in 
the network. The subscript A mean that the correspond- 
ing heuristics is applied to the nodes from the commu- 
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nity A. In the results presented below we used commu- 
nities of size Na = 1000, Nt = 4000, and p^n = 0.8, 
and the degree distribution is characterized by 7 — 2.5, 
and the lower cutoff is m = 10. Each point was aver- 
aged for 100 random trials. One can see that for small 
M, the difference between two targeting strategies are in- 
deed significant, similar to the results in random Erdos- 
Renyi graphs. Note that in absolute terms, the difference 
for the maximum degree heuristics is significantly larger 
compared to the greedy targeting case. This difference 
however, diminishes as one increases the cutoff. Thus, for 
networks with very large cutoffs, targeting nodes from the 
smaller group does not provide any improvement. In fact, 
depending on the actual cutoff mechanism, one might 
be better of targeting nodes from the larger community. 
For instance, for networks with a power-law exponent 7, 
the so called natural connectivity cutoff scales with the 
network size as ^ 

Ari/(7-i) [2]: Thus, one might argue 
that larger community should have a higher cutoff, which 
might make it more beneficial to target nodes from that 
community. 



IV. DISCUSSION 



mogeneous structure, then a newly introduced idea can- 
not survive unless it is initially adopted by a large num- 
ber of individuals. However, if the network describing the 
interaction of sub-units meets certain structural condi- 
tions, then the idea might take over the whole population 
even starting from a small number of initial adopters. 
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APPENDIX 

We first focus on the random targeting case. Let us 
define S — -f — 2, and (3 = m/M, and x = mO. Then the 
steady state Equation |8] corresponds to the zeros of the 
following function: 
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In this paper we examined the problem of maximizing 
influence propagation in structured heterogeneous net- 
works. We demonstrated that for models with critical 
behavior, the structural properties of the network, and 
specifically, its community structure, can have impor- 
tant implications for the influence maximization problem. 
We demonstrated analytically that for two-community 
networks, targeting nodes from the smaller community 
might cause a global cascade with significantly fewer 
number of seed nodes. This effect becomes increasingly 
important if the sizes of two communities are vastly dif- 
ferent. We also showed through numerical simulations 
that a similar picture holds for multi-community net- 
works. 

In practice, one of course does not have precise esti- 
mates of model parameters, such as activations thresh- 
olds, or even the precise topology of the network. Thus, 
the problem of finding the optimal target set for influence 
maximization might not be well defined. On the other 
hand, with the surge in online networking sites, infor- 
mation about underlaying community structure in such 
networks is often available. Our results suggest that in 
such cases, paying attention to the community structure 
of the network might be beneficial for influence maxi- 
mization. 

We also note that the networks considered here mimic 
scenarios where innovations are introduced through a 
small community of early adopters. In this respect, our 
work resonates well with organizational viscosity model 
of Krackhardt [211 US] th&t describes diffusion of ideas in 
an organization. Here organization is modeled as a num- 
ber of interacting sub-units, with closer social ties within 
each unit. When the organization has a more or less ho- 



A simple inspection shows that x = to is always a solution 
(aside from exponentially small corrections). For 5 > \ 
(7 > 3), there is a critical fraction pc below which other 
solution appears, as it is schematically shown in Figure [4] 
Thus, wc need to show that for (5 > 1, Equation |A. 1| has 
a solution for small x. Note that derivative of g{x) is 
given by 



(A.2) 

Consider the case /? = (M 00). Starting at g{x — 
0) ~ po, the function g{x) is a strictly decreasing over 
the interval < x < xq, where xq is determined from 
g'{xo) = 0, which yields 
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There are two separate cases: For H—5 > 0, the integral, 
for small xq, can be replaced by the gamma function 
T{H — 5), which yields 



{H-iy. 
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And for H — 5 < 0, after integrating by parts and keeping 
the leading term we obtain 



{6-H){H-iy. 1 



1 - Po 



(A.5) 



In both cases, xq remains finite as po — > 0. Thus, for 
sufficiently small x, one has g{x) w pg — cx, where c > 
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does not depend on po- Consequently, g{x) will intersect 
zero for sufficiently small po- More precisely, the critical 
fraction pq = for which the other solution appears is 
found from 5(2:0) — 0. 

Now consider the case < (5 < 1, (or 2 < 7 < 3). 



Without connectivity cutoff {(3 = 0) the integral in A. 2 
X for small x does not depend on x, and can be approx- 
imated by V{H — 5). Thus, the derivative g'{x) behaves 
as g'{x) ^ \/x^~^ for small x. Consequently, there is 
no other solution except for the one at a; = m. A similar 
argument holds for 7 = 3, where g'{x) remains finite, but 
positive for small x. The situation changes as one intro- 
duces a finite cutoff M . Indeed, for finite, but small /3, 
/3 <C 1, the integral is approximately {x/(3)^~^ /{H — 5). 
Thus, the derivative is negative over a finite interval 
< a; < xq, where 
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According to the same argument as above, there is a 
Po = Pc such that g{xo) — 0. 

Let us now consider the maximum degree heuristics, 
for which we have 



g{x) = 1 - - - -^x' / dte-H-'-^ V 



Here (3^) = m/ ttiq , and the cutoff connectivity too is given 



by Equation [9] which, for f3 <^1, reads 

r, _ 1/(1 + 5) 

Po — Po 
It is easy to check that 
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Furthermore, the derivate of g(x) is 
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Consequently, g(x) is negative for < x < Xq where 



H-S 



iS-H){H-l)\ 



Thus, when decreasing po, the interval where 17(0;) de- 
creases shrinks as with ai = -jjjz^jjj^rg^ ■ At the 
same time, g{x = 0) = p^^, with 012 = 5 /{I + 6). For 
6 > 1 one has ai > a2, which suggests that g{x) will 
cross the zero at some critical value po = Pc- And in con- 
trary, for (5 < 1 one has ai < a2, which means that g(x) 
always remains positive as po — ^ 0; thus suppressing crit- 
ical behavior. Finally, repeating the arguments above, 
one can show that introducing a connectivity cutoff for 
6 < Iwill recover the criticality . 
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